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A Josephson junction may be driven through a transition where the superconducting condensate 
favors an odd over an even number of electrons. At this switch in the ground-state fermion parity, 
an Andreev bound state crosses through the Fermi level, producing a zero-mode that can be probed 
by a point contact to a grounded metal. We calculate the time-dependent charge transfer between 
superconductor and metal for a linear sweep through the transition. One single quasiparticle is 
exchanged with charge Q depending on the coupling energies 71,72 of the metal to the Majorana 
operators of the zero-mode. For a single-channel point contact, Q equals the electron charge e in 
the adiabatic limit of slow driving, while in the opposite quenched limit Q = 26 .^ 7172/(71 J- 72 ) 
varies between 0 and e. This provides a method to produce single charge-neutral quasiparticles on 
demand. 


I. INTRODUCTION 

Superconductors connected by a Josephson junction 
can freely exchange pairs of electrons, but single-electron 
transfer is suppressed by the superconducting gap.® The 
tunneling of an unpaired electron into the junction is 
an incoherent, stochastic source of charge noise in a 
Cooper pair transistor.® In contrast to this undesirable 
“quasiparticle poisoning”, a controlled phase-coherent 
way to exchange single quasiparticles with a supercon¬ 
ductor would be a desirable tool, that would comple¬ 
ment existing single-electron sources in normal metals 
and semiconductors.®® 

Here we propose to exploit the phenomenon of a 
fermion-parity switch to transfer, phase coherently and 
on demand, a single quasiparticle of adjustable charge 
Q from a Josephson junction to a metal probe (see Fig. 
&)■ A fermion-parity switch is a topological phase tran¬ 
sition (zero- dim ensional class D in the “ten-fold way” 
classificatiorPEni) where the superconducting condensate 
can lower its ground-state energy by incorporating an un¬ 



phase difference <p 


FIG. 1: (a) Josephson junction formed by a superconducting 
ring interrupted by a nanowire. The junction contains two 
Majorana zero-modes, separated by a tunnel barrier (height 
Vo). A time-dependent flux 'l>(t) through the ring drives the 
phase = 4>(t) X 2e/?i through a fermion-parity switch, 
at which a quasiparticle is injected as a current L{t) into the 
grounded metal probe, (b) Pair of phase-dependent Andreev 
levels ±Eo{(f>) in the closed Josephson junction (uncoupled 
from the metal). The switch in the ground-state fermion par¬ 
ity vf is signaled by a level crossing. 


paired electron and changing the number of electrons in 
the ground state from izp even to zzp odd,^ leaving behind 
as “defects” an odd number of quasiparticle excitations 
above the ground state. 

In the quasiparticle excitation spectrum, the switch 
in the ground-state fermion parity is signaled by the 
crossing of a pair of bound states (Andreev levels) at 
E = 0 (the Fermi level). There may be an even num¬ 
ber of switches when the phase difference (f) across the 
Josephson junction is incremented by 27r — if there is an 
odd number of switches (as in Fig. the superconduc¬ 
tor is topologically nontrivial. The two lowest Andreev 
levels ±Eo{(j)) of a nontrivial Josephson junction have a 
cos((/)/2) phase dependence, 

Eoi(f) = Aoy^oCOsW2). ( 1 ) 

The superconducting gap is Aq and Tq G (0,1) is the 
transmission probability through the junction. For small 
To this describes a pair of bound states at nearly zero 
energy, consisting of an equal-weight superposition of 
electron and hole excitations. Such a charge-neutral 
quasiparticle is called a “Majorana fermion” (or Majo¬ 
rana zero-mode) because of the identity of particle and 
antiparticle. These objects have unusual non-Abelian 
statistics (see Refs. mm for recent reviews), but here 
it is only their charge-neutrality that matters. 

Fermion-pariW switches are actively studied, 

theoreticalljff^^^ and experiment ally, I 22 H 2 II for the 

connection to topological superconductivity and Ma¬ 
jorana fermions .12SH2SI The dynamics of the transition 
is what concerns us here, in particular the quench 
dynamics, where 4i{t) is driven rapidly through the 
switch from even to odd ground-state fermion parity. 

The geometry of Fig. that we consider is modeled 
after existing experiments {e.g., Ref. HI]), where a meso¬ 
scopic Josephson junction is formed by a semiconduc¬ 
tor nanowire connecting two arms of a superconduct¬ 
ing ring. A time-dependent flux $(t) enclosed by the 
ring imposes a time dependence on the phase difference 
(j){t) — $(t) X 2e/Ii across the junction. When the Joseph¬ 
son junction is quenched through a fermion-parity switch 
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there will appear a current pulse I(t) from the supercon¬ 
ductor (S) into the metal (N). We seek the quasiparti¬ 
cle content of that pulse. How many quasiparticles are 
transferred? What is the transferred charge? In partic¬ 
ular, we wish to establish the conditions under which a 
single quasiparticle is transferred with vanishing charge 
expectation value. 

We find that the quench dynamics transfers one single 
quasiparticle from the superconductor to the metal, as a 
wave packet that is a coherent superposition of electron 
and hole states near the Fermi level. A nearly charge- 
neutral equal-weight superposition is produced in a topo¬ 
logically nontrivial superconductor, if the metal probe 
couples predominantly to one of the two spatially sepa¬ 
rated Majorana zero-modes. More generally, for two ar¬ 
bitrary coupling constants 71,72 we derive that the quan¬ 
tum quench injects a charge 

Qquench — 2cy^7lb^/(7l A 72 ) (2) 

into a single-channel point contact. For a multi-channel 
point contact the injected charge is reduced further by 
a factor TZ determined by the peak height Gpeak = 
(4e^//i)(l — TZ^) of the point contact conductance at res¬ 
onance. 


II. MICROSCOPIC MODEL 


Before proceeding to the mathematical analysis of the 
quench dynamics, we explore the relevant physical pa¬ 
rameters in a microscopic modeP^ for an InSb nanowire 
(length L = 2.5 ^m, width W = 0.25 ^m, Fermi energy 
Ep = 1.52 meV, corresponding to 4 occupied electron 
subbands), coupled at both ends to a Nb superconductor 
(induced gap Aq = 0.4 meV). Spin-rotation symmetry 
is broken by Rashba spin-orbit coupling (characteristic 
length 4 o = ImoSCtso = 0.25p,m), and time-reversal 
symmetry is broken by a magnetic field parallel to the 
wire (Zeeman energy Vz = ^geSiJ-BB = 0.6meV). For 
these parameters, the Josephson junction is in the non¬ 
trivial regim e, wit h a pair of Majorana zero-modes at 
the two ends.l^i^ We tune the coupling strength of the 
Majoranas by means of a tunnel barrier of width 25 nm 
and adjustable height Vq (which might be experimentally 
realized by means of a gate voltage). The data shown in 
Fig. i is for Vo = 15meV. (See App.|^for details of the 
calculation.) 

The Josephson junction is coupled by a point contact 
to a normal-metal probe, which plays the role of a fermion 
bath that can exchange quasiparticles with the super¬ 
conductor. We assume that the charging energy of the 
junction is much smaller than the Josephson energy, to 
ensure that the Coulomb blockade of charge transfer is 
not effective. The Josephson junction is now an open 
system, with quasibound Andreev states — ir„ that 
acquire a finite life time h/2Tn. The evolution of a pair of 
these states through the fermion-parity switch is shown 



FIG. 2: Phase dependence of the complex energies En — iVn 
of a pair of quasibound states of the open Josephson junction 
(solid curves), when the energies ±Eo of the closed junction 
(dashed curves) vary through the level crossing of Fig. [^. 
At the fermion parity switch, the inverse lifetimes r„ reach 
opposite extremal points TV'yn, n = 1,2. 

in Fig. [2p^ The coupling constants 7 „ that determine the 
transferred charge can be read off from 

TTJri = lim r„(0). (3) 

(f)—¥TZ 

Particle-hole symmetry requires that the complex en¬ 
ergies come in pairs ±E — iF, symmetrically arranged 
around the imaginary axis. This constraint produces 
a bifurcation point (pole transitiorl^ or exceptional 
poinlP^ at which the real part is pinned to i? = 0 and the 
decay rates Fi, r 2 become distinct — resulting in widely 
different 71 , 72 . The unusual extension of the level cross¬ 
ing over a finite interval seen in Fig. [^is the key distin¬ 
guishing feature of level crossings in superconducting and 
non-superconducting systems, and makes the dynamical 
problem considered here qualitatively different from the 
familiar Landau-Zener dynamics.l^ 


III. SCATTERING FORMULATION 

The exchange of quasiparticles across the NS interface 
is described by the scattering matrix 

S{t, t') = 5{t - t’) - 27rfW^G(t, t')W. (4) 

The coupling matrix W to the fermion bath is assumed 
to be time-independent. The retarded Green’s function 
G{t, t') satisfies the differential equatioiP^ 

{id/dt-H[(j){t)]+i-KWW'^)G{t,t') = 5{t-t'), (5) 

where E[{(j)) is the Bogoliubov-De Gennes Hamiltonian of 
the Josephson junction at a fixed value (j) of the super¬ 
conducting phase difference. (We have set h = 1 for ease 
of notation.) Fourier transform to the energy domain is 
defined by 

/ OO pOO 

dt / dt' ( 6 ) 

-OO J —00 
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In a stationary situation, with a time-independent 
Hamiltonian H, the scattering matrix is diagonal in en¬ 
ergy, S{E,E') = 2'k 5{E — E') So{l^ with Sq given by 
the Mahaux-Weidenmiiller formula,*^ 

SoiE) = 1 - 2TnW^(E - H,s)-^W, 
HeS=H-i-KWWK 


The formulation of this dynamical problem in an open 
system in terms of an effective non-Hermitian Hamilto¬ 
nian HeS go es back to the early days of nuclear scattering 
theoryP^® 

For a minimal description, we take a pair of An¬ 
dreev levels in the Josephson junction coupled to a pair 
of electron-hole modes in a single-channel metal probe. 
(The multi-channel case is addressed in Sec. VB| ) Both 
E[ and W are now 2x2 matrices. Particle-hole symmetry 
requires that 


H = W = ( 8 ) 


(The Pauli matrix interchanges electron and hole in¬ 
dices.) Particle-hole symmetry is the only symmetry con¬ 
straint we impose on the system (symmetry class D), 
assuming that time-reversal symmetry and spin-rotation 
symmetry are both broken by magnetic field and spin- 
orbit coupling in the nanowire. 

Using also that iJ = we have the general form 


H = Eoa^, IT = 



with real coefficients a, a', X±. The eigenvalues 71,72 > 
0 of the coupling matrix product WW^ are given by 

7 i = (^+ + ^-)^) 72 = (A+-A_)^. (10) 

The eigenvalues of iJeS (representing the poles of Sq 
in the complex energy plane) are given by 

E± = -inj ± Eo\/l + {tt^/EqY - {nj/Eo)^, ( 11 ) 

in terms of the arithmetic and geometric mean 

7 =^ 71 + 72 ), 7 = V7i72- (12) 

The evolution of E± through the fermion-parity switch 
is shown in Fig. The relation E+ = —EZ. required by 
particle-hole symmetry produces a bifurcation point at 
whic h the two quasibound states acquire distinct decay 
rates,see also Fig. 

The time dependent phase difference across the 
Josephson junction shakes up the fermion bath in the 
normal metal. We assume zero temperature, so that the 
unperturbed Fermi sea is the vacuum state |0) for exci¬ 
tations: a(if)| 0 ) = 0 for if > 0 , with a = ( 01 , 02 ) the 
two-component Nambu spinor of annihilation operators 
for Bogoliubov quasiparticles. The fermion-parity switch 
produces a superposition 


(13) 




1 


1 —Try -|- Try/y^ — y^ 

1 E_ 

E+ 1 


—Try — 7r\/7^ — y^ 1 


0 ReE± 


FIG. 3: Evolution of the complex eigenvalues E± of the effec¬ 
tive Hamiltonian 0 of the open Josephson junction (coupled 
to a metal probe), when the real eigenvalues ±i?o of the closed 
junction vary through a level crossing. At the fermion parity 
switch, if 0 = 0 and E± reach opposite extremal points on the 
imaginary axis. 


of the vacuum state with p-particle excited states 




^ ^ aHE)S{E,E')aiE') 


.E>0E'<0 


|0). (14) 


(The sum J2e evaluated as (27r)“^ / dE.) The weight 
(0 of the unperturbed Fermi sea follows from the normal¬ 
ization (ff'l'k) = 1 . 


IV. LINEAR SWEEP THROUGH THE 
FERMION-PARITY SWITCH 


We now proceed to a complete solution of the dynam¬ 
ics of the fermion-parity switch, to derive the result 0 
for the charge of the transferred quasiparticle. The non¬ 
superconducting counterpart to this problem was studied 
by Keeling, Shytov, and Levitov.SH Their analysis pro¬ 
vided much guidance and inspiration for what follows. 

We calculate the scattering matrix for a linear sweep 
through the fermion parity switch: EQ[(j>{t)] = Re¬ 
ferring to Eq. 0. this linear approximation of the spec¬ 
trum is justified for rapidities 7 q ^ VTo^oy- In the 
energy domain, Eqs. Q and ([^ then take the form 

S{E,E') = 2tt6{E-E') - 27rie-*“'""AG(E,E;')Ae*“'"7 
{i-fla^d/dE + E + mK^)G{E, E') = 2 'k5{E - E'). 

(15) 


The solution for the Green’s function factorizes, 


G{E,E’) = 
Q{E - E') 


^X{E)&{E - E')a,X-\E')a,, 
*7o 

_ (6{E-E') 0 \ 

“ 0 e{E’ -E)) ■ 


(16) 


l'h) = Co|o) + Er=il'^'p) 


(17) 




















4 


Here 0{E) is the unit step function and the matrix X{E) 
solves the homogeneous equatioiP^ 

{i-ila^d/dE + E + mK^)X{E)=Q. (18) 

Because of particle-hole symmetry, X has two rather 
than four independent elements, 


X{E) = 


fu{E) v*{-E)\ 
[v{E) u%-E)) ’ 


(19) 


determined by 

7 qm" -I- (e^ + — i)u = 0, 8v = ieu — 7oM^ (20) 
e = (£; -f t7r7)/7o, 5 = - ^ 2 )ho- ( 21 ) 


The retarded Green’s function is specihed by G —?► 0 
in the limits E —?> -l-oo or E —>■ — 00 . The factor 0 in 
Eq. (16 1 ensures that this two-sided decay follows from 
the one-sided decay u, z; —>■ 0 for E —-|-oo. With this 
condition the solution of Eq. (20) read^^ 


u{E) = U 

v{E) = e*® — 5 ; —ie^), 


( 22 ) 


where U is t he co nfluent hypergeometric function of the 
second kind.l^^*^ The determinant of X is particularly 
simple (see App. H 

DetX = exp(— 7 r^^/ 4 ), (23) 



7/7 


FIG. 4: Expectation value of the charge of the quasiparticle 
transferred between the superconductor and a single-channel 
metal probe, following a fermion-parity switch with rapidity 
70 . The charge Q is given as a function of the ratio 7/7 
of the geometric and arithmetic mean of the coupling ener¬ 
gies to the two Majorana operators involved in the transition. 
The curves are calculated numerically from E q. (|28| . The 
quenched and adiabatic limits are given by Eqs. (|29[l and (|32[|. 


with wave function z/i given by Eq. (25). The transfer 
of this quasiparticle is observable as an electrical current 
pulse, with expectation value 

= (27) 


independent of energy. 

The scattering matrix (15) results as the dyadic prod¬ 
uct of two vectors, 


S„miE, E')\e>E' = -hiEWmi-E'), (24) 

hE) = (2^/7o)e"^^/«e—. (25) 


Substitution into Eq. (14) gives Idfp) = 0 for p > 2 be¬ 
cause of the anticommutation of the creation o per ators, 
so that only a single-particle excitation remains,S£l 


The expectation value of the total transferred charge Q = 
I{t)dt is given by 

Q = ^(A^ - A2_) r dE {\u{Eh - HEh). 
7o Jo 

(28) 

For definiteness we take A^ > A?, in what follows (other¬ 
wise the sign of currents and charges should be inverted). 


V. TRANSFERRED CHARGE 


Ivki) = - 5] ^ [hE)aHE)][ri-E')aiE')m- (26) 

E>0 E'<0 

This absence of multi-particle excita tions is a generic fea¬ 
ture of rank-one scattering matrices 

The normalization X]£;>o IV’(A)P = 1 can be derived 
directly from Eq. ( |I^ . (See App. i) This implies that 
(fkijlki) = 1 , hence there is no contribution from the 
vacuum state [(Co = 0 in Eq. Q]- Corrections of order 
I = exp(—27ri?7/7Q) to the normalization appear be¬ 
cause of the finite band width E < -^/To Ag. Since we have 
assumed 7 o ^ t/Tq A 07 we can ascertain that the sweep 
through the fermion-parity switch will fail to produce a 
quasiparticle with exponentially small probability. 

The Josephson junction thus injects a single Bogoli- 
ubov quasiparticle into the metal probe, in a pure state 


A. Single-channel probe 


A single quasiparticle passes through the NS interface 
irrespective of the rapidity 70 , but the transferred charge 
differs. Fig. shows results from a numerical evaluation 
of Eq. (28). Analytical results can be obtained in the 


quenched limit 70 7 i, 72 of a fast fermion-parity switch 
and in the opposite adiabatic limit 70 ^ 7 i, 72 of a slow 
switch. 

In the quenched limit we set 5 —> 0 and since 
[7(0,1; —ie^) = 1 we have u —>■ exp(i£^/2), u —)■ 0. The 
current and transferred charge evaluate to 

7quench(7) = 27re7 exp(-27r7t)6»(t), Qquench = 07 / 7 . 

(29) 

This is the result ([^ announced in the introduction. 
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The adiabatic limit may be obtained, with some ef¬ 


fort, from the Fourier transform (27) in saddle-point ap¬ 


proximation, or more easi ly by starting directly from the 
general scattering formulcPSHSI] 


4diabatic(t) = ^ TrS'|,(0,t)CT2—S'F(0,t). 


47r 


(30) 


(A selfcontained derivation of this formula is given in 
App. [^) The adiabatic charge transfer is described by 
the “frozen” scattering matrix 


5F(A,t) = 5o(A)|^^^p), (31) 

with So from Eq. ([^ evaluated for a fixed value of 
the phase across the Josephson junction. The result is 


-^adiabatic (t) — 


V7i72 


TT^yi72/70 + 7o^^ ’ 


^adiabatic — C. (32) 


The exponential versus Lorentzian current prohles 


(29) and (32) have the same form as in the non¬ 


superconducting problem of Ref. 1411 but there the trans¬ 
ferred quasiparticle was an electron of charge e. Here 
what is transferred is a Bogoliubov quasiparticle, which 
is not in an eigenstate of charge. In the quenched limit Q 
can vary between 0 and e, depending on the ratio of the 
geometric and arithmetic mean of the two coupling ener¬ 
gies 7 i, 72 of the metal probe to the Majorana operators 
of the zero-mode. A nearly charge-neutral quasiparticle 
is transferred if 71 <C 72 , when Q = 26 -^/ 71/72 in the 
quenched limit. 


B. Multi-channel probe 

So far we have assumed that the metal probe supports 
a single electron-hole channel. More generally, the cou¬ 
pling between the superconductor and the metal would 
involve N electron-hole channels, where N would include 
both orbital and spin degrees of freedom. This multi¬ 
channel generalization is worked out in App. A single 
quasiparticle is injected, as before, with a reduced charge 
Qn = Ti-Qi- The reduction factor TZ G [0,1] is indepen¬ 
dent of the rapidity 70 . It is determined entirely by the 
point contact conductance, which at the fermion parity 
switch has a resonant peak of height 

Gpeak= ^(l-7^2)■ (33) 


a charge expectation value that can be adjusted between 
0 and e. A nearly charge-neutral quasiparticle is pro¬ 
duced in the quenched limit of a fast parity switch, if 
the metal couples predominantly to a single Majorana 
operator in the Josephson junction. This device could 
be used for superconducting analogues of single-electron 
collision experiments,^^® such as the Hanbury-Brown- 
Twiss or Hon g-Ou-Mandel interferometer for Majorana 
fermions 

Experimentally, one can determine the value of Q by 
sweeping up and down through the fermion-parity switch 
and measuring the shot noise power Pshot- In each pe¬ 
riod T a charge { 0 , -|-e, —e} is transferred with probability 
{1 - 2p(l -p),p(l -p),p(l -p)}, where Q/e = |1 - 2p| 
is the average charge transferred during a sweep up or 
down. The full distribution of the transferred charge is 
trinomial. The first moment vanishes and the second 
moment is given by 

Pshot = 2p(l -p)(e^/r) = |r“^(e^ - Q^). (34) 

Referring to the model calculation of Eig. a band 
width of VTqAo — 10 GHz at a driving frequency of 
1/r ~ 0.1 GHz would imply a rapidity 70 — 1 GHz (so 
that 'JqT ~ \/ToAo)- The escape rate 7 could then vary 
between, say, 0.2 GHz and 2 GHz to vary between the 
adiabatic and the quenched regime. These frequencies 
should all lie above the decoherence rate of the Bogoli¬ 
ubov quasiparticle due to charge noise, which could be 
below 1 MHz.l® 

An alternative way to measure the transferred charge 
is to apply a voltage V between the two superconductors. 
The phase will then advance with constant rate d(f>/dt = 
2eV/h, producing a current I = Q x 2eV/h (assuming a 
single level crossing in a 27r phase interval). 
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Appendix A: Model Hamiltonian 


VI. CONCLUSION 


The model Hamiltonian for the nanowire Josephson 
junction of Eig. has the Bogoliubov-De Gennes form 


In conclusion, we have investigated the phase-coherent, 
deterministic counterpart of incoherent, stochastic quasi¬ 
particle poisoning: A fermion-parity switch in a Joseph¬ 
son junction transfers a single quasiparticle into a metal 
contact, on demand and in a pure state. The quasiparti¬ 
cle is a coherent superposition of electron and hole, with 



-avH^{-p)ay) ’ 


(Ala) 


Ho = 


2mes 


Ef + -^{(JxPy - (TyPx) + \gefiPBBax 


+ Vo [0(x - VFb/2) - 0(1 - IFb/2)]. (Alb) 
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Gpeak = 3.97eV/i 


An = 0.4 meV 


E-p = 1.52 meV e®‘^An 
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B 


Vz = 0.6 meV 


FIG. 5: Nanowire Josephson junction modeled by the Hamil¬ 
tonian (All, discretized on a square lattice (lattice constant 
a = 25 nm). The InSb nanowire is grey, with a tunnel barrier 
(width 25 nm) in black, the superconducting contacts are yel¬ 
low, the normal-metal probe (width 100 nm) is blue. There 
are 4 electron subbands in the nanowire and 8 in the probe, 
counting spin. The peak conductance at the fermion-parity 
switch is indicated. 


Electrons and holes are coupled by the induced s-wave 
pair potential A at the superconducting contacts, with 
a phase difference (j). The single-particle Hamiltonian 
Hq contains Rashba spin-orbit coupling and the Zeeman 
energy of a magnetic field parallel to the nanowire. A 
potential barrier of strength Vq and width Wb is located 
at the center of the junction. 

The Hamiltonian H is discretized on a square lattice, 
to obtain a tight-binding model.l^For the parameters in¬ 
dicated in the figur e, the Josephson junction is in the non¬ 
trivial regime,with a pair of Majorana zero-modes at 
the normal-superconducting (NS) interface, weakly cou¬ 
pled via the potential barrier. A normal-metal lead is 
attached perpendicular to the nanowire, coupling pre¬ 
dominantly to one of the two zero-modes. 

To obtain the complex energies of the quasibound 
states, the imaginary part of the lead self-energy is added 
to the tight-binding Hamiltonian of the junction. Diago- 
nalization of this non-Hermitian Hamiltonian yields the 
complex eigenvalues En[4>) ~ *r„(^) plotted in Fig. 


Appendix B: Details of the calculation of the 
Green’s function 


u, V solve 


{i^lazd/dE + E + inA^) 


= 0 . 


(B2) 


This allows us to express the derivatives 


vS) ' 


u(E)\ 

v{E)J ’ 


(B3) 


_ ^i-E) 

dE \v*{-E) 


= - yaz{E + 

7o 


u*{-E)\ 
v*{-E) ■ 


(B4) 


Since A is a real and symmetric matrix, it follows that 

T 


d i 

— Det A = ^ 

dE % 


u*{-E) 

v*{-E) 


[{E -k mA^) 




= 0 , 


(B5) 


so Det X is independent of E. 

From Eq. (22 1 we have an explicit expression for the 
determinant of X : 

Det A = —ie^)U 


- i; + \^5\ Eie^). (B6) 


This is an analytic function of e = (F-l-iTryj/yo, which is 
independent of E and hence independent of e. At £ = 0 
we may evaluate it by means of the identitie^l 


C/(a, i,0) = 


r(i+a)’ 


r(i + la)r(i — ia) = 

r(l + ia)r(l — ia) = 


cosh na ’ 

TTO 

sinh 7ra 


(B7) 


(B8) 


Substitution into Eq. (B 61 at e = 0 gives 


1. Evaluation of the determinant 


Det A = exp(— 7 ri 5 ^/ 4 ), 


(B9) 


Since the expression (ITm for the Green’s function con¬ 
tains both the matrix A^F) and its inverse, we need to 
evaluate the determinant of this 2 x 2 matrix. As a first 
step we will show that Det A is energy independent. This 
can be done directly from the differential equation (18) 
for A. 

We write the determinant in the form 


Det A(£;) = 


u*{-E) 

v*i-E) 


u{E) 

v{E) 


(HI) 


as in Eq. (231. 


2. Normalization of the excited state 


We wish to demonstrate that the wave function (25) 
of the single-particle excited state is normalized to unity. 
For that purpose we need to evaluate the integral 




2'KdE fu*{E) 


lo -f^BetX \v*iE) 


A^ 


<E) 

v{E) 


(BIO) 


and take the derivative with respect to E. The functions 
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We again use the fa ct that u, v solve Eq. ( |B2[ ) . Substi¬ 
tution into Eq. (BIOI gives (denoting u' = du/dE) 


I{t) = (lit)), given by 


A7 = 


-2 

DetX 


POO 

/ dE [u*u'— v*v'— iEjQ‘^(uu* + vv*)] 

Jo 


+ 


5 ^(l»( 0 )P-|..( 0 )| 
2 


DetX 


= 2-AT* 


pOO 

/ dE [uu*' — vv*' + iEjQ^{uu* -|- vv*)] 

Jo 

(Bll) 


and because Af is real, we indeed have AC = 1. Notice 
that = 1 also implies = 1 in Eq. (26). 


Appendix C: Scattering formula for the charge 
transfer in the adiabatic regime 

The current passing through the NS interface in the 
adiabatic regime 70 ^ 71 , 72 of a slow fermion-parity 
switch can be evaluated most easily from the scatter¬ 
ing formula ( |30| , which is the analogue for Bogoli- 
ubov quasiparticles of a well-known formula for normal 
electrons For completeness we give a derivation of 
Eq. (|3§. 

One subtlety in this derivation is that Fourier trans¬ 
forms of quasiparticle annihilation operators a{E) to the 
time domain need to include both positive and negative 
energies in order to produce a complete basis set. This 
results in a double counting of the quasiparticle excita¬ 
tions, because of the relation a{—E) = axa'^{E). To cor¬ 
rect for the double counting we include a factor 1/2 in 
the definition of the current operator,!^ 


I{t) = 

aontit) = / — e ^^aont{E). 

7-00 


(Cl) 


The outgoing and incoming operators are related by 
the scattering matrix, 

%-S{E,E')aUE'), (C2) 

-00 

which satisfies the unitarity condition 

— Y,Snn'{EuE')S*^^,{E2,E') 

- 2^ (C3) 

= 27r5„m(5(Ei - E 2 ). 

The incoming operators have the equilibrium expecta¬ 
tion value 

{aiiE)am{E')) = 2tt 5{E - E')Snmf{E), (C4) 

with f{E) = (1 -I- the Fermi function at tem¬ 

perature T. We seek the current expectation value 


m)=^-e r r — r 

2 27r 2n J 27r 

X f{E')TiS\E + uj,E')a^S{E,E'). (C5) 


Because of the unitarity condition (C3), the integral 


over E' without the factor f{E') vanishes, 
dE' 


27r 


Tr S'1' {E + w, E')a^S{E, E') = 2tt5{uj) Tr ( 


= 0 . 


(C 6 ) 


We may therefore equivalently write 

m) = ie r d_E r r 

2 2 ^ 27r 27T 

X [fiE') - f{E)] Tr (ij + E')a,S{E, E'). 

(C7) 

It is convenient to introduce the Wigner transform 

— e-'^'^S{E+\E',E-\E'), (C 8 ) 

because it becomes the frozen scattering matrix Sp{E, t) 
from Eq. © in the adiabatic limit .1^ More precisely, 

Sw{E+SE,t) = SFiE,t)+Oi-fo/E,)+0{6E/E,), (C9) 

with Ec = min( 7 i, 72 ) the width of the quasi-bound 
state. 

Fourier transformation of the time variable gives 


Svj{E,uj)= / dtE‘^*Sw{E,t) = S{E + ^uj,E - ^uj). 

J —00 

(CIO) 

In terms of Sw{E,uj) the expression (C7) for the current 
reads 


,, , , f°° dE duj' 


duj 

27r 


X [f{E-\u:')-f{E+\u:')] 

X Tr5'^(i/ -I- 50 ;, w -I- uj')azSw{E,oj'), (Cll) 

with the definitions E = ^{E + E'), uj' = E — E'. 

The integrals over uj and to' contribute over the range 
~ 7 o uj,uj' < 7 o. To leading order in 70 we therefore 
have 

Tr S'^(i? -I- Iw, oj + uj')azSw{E, oj') = 

Tr4(E,w + w')(T,5F(E,w')+t7(7o/^c), (C12) 


in view of Eq. (|C9|). Substitution into Eq. (Cll|), with a 












change of variables oj” = oj + ui', results in 


dE 


du;' 


duj'‘ 


m = he ^ ^ 




27 T 


27T 


27 T 


X [/(^-y)-/(^+y)] 

X TrSl{E,u;")a,SFiE,uj')[l + O{jo/E,)] 


= b 


■'F\ 

dE 


2tt 


dio 

27r 


X [f{E-lu;)-f{E+luj)] 

X Tt sUE,t)a,SFiE,io)[l + O{jo/Ec)]. (C13) 

Since we do not wish to assume that 70 is small com¬ 
pared to kT, we expand the difference of Fermi functions 
in square brackets to all order in w, 

^ (2p+l)! dE^P^ ^ ^ 

I / y I 1M a7:i9r> Oj.2r) I J (EjUje 


\p=0 


{ 2 p+iy.dE‘^P dt^P 


(C14) 


Upon partial integration, the sum over p contributes to 
the integral (|C13|) terms of order 


g2p Q2p 

O&pWp 


SF{E,t)=0{jo/EcfP, 


(CIS) 


of electron-hole modes in the normal-metal probe. This 
coupling is described by the 2x2 coupling matrix W 
defined in Eq. ©• More generally, a multi-channel probe 
has a 2 X 2N coupling matrix of the form 

w = (Wi,W 2 ,... Wn), W„ = (^2 ai) ’ 

constrained by particle-hole symmetry: W = axW*ax- 
We collect the complex coefficients an, Pn in a pair of 
vectors. 


a = (ai, 02 ,... ow), P = {Pi,P 2 , ■ ■ ■ Pn), 
and define the inner products 

N N 

(a|a) = ^ |a„|2, (/3|/3) = ^ |/3„p, 

n—1 n—1 

N 

(«l/3) = '^alPn- 

n—1 

The decay rates 71 , 72 of the pair of quasibound An¬ 
dreev levels are the eigenvalues of the 2 x 2 matrix 

N 

WW^ = 

n—1 


(D2) 


(D3) 


f(ala) + (m 2(a|/3)* \ 

2(a|/3) (ajaj + imj’ 


(D4) 


71 = (ala) -h (/3|/3) -P2|(a|/3)|, 

72 = (a|a) -h (/3|/3) - 2|(a|/3)|. 


(D5) 


so only the p = 0 term needs to be retained to leading 
order. 

We thus arrive at 


duj 


r ^ r 

2 27T 

X Ti SliE,t)a,SFiE,io)[l + Oi-fo/Ec)] 


nt)=-hl ^/'(E)a.e—‘ 


r°° dE 


271 


xTTSl{E,t)a,-^^SFiE,t)[l + 0{jo/E,)]. 


(C16) 


At zero temperature, when —f'{E) 5{E), we recover 

Eq. (pOl), 


^e „ 9 


4diabatic(t) = — Tr 51,(0, t)cr^—S'F(0,t). 
Att at 


(C17) 


Appendix D: Multi-channel probe 
1. Coupling matrix 


As before, we define the arithmetic and geometric aver¬ 
ages, 

7 = 1 ( 71 + 72 ), 7 = V7i72- (D6) 

For later use, we also note that 

N 

Wa,W^ = = ((a|a) - (/3 |^))(t,. (D7) 

n—1 


2 . Scattering matrix 


Carrying through the same steps as in the single¬ 
channel case, we have the following expression for the 
27V X 27V scattering matrix S in terms of the 2x2 Green’s 
function G: 


S{E, E') = 2tt 5{E - E') - 2'KiW'^G{E, E')W, 

[ 1 - 11 ( 7 ^-^+E + mWW"^^ G{E,E') = 2 ttS{E-E'). 

(D 8 ) 


The solution for G has the factorized form (16), in terms 
of the 2 x 2 matrix 


In the main text we assumed that the pair of Andreev 
levels near the level crossing is coupled to a single pair 


X{E) = 


(u{E) v*{-E)\ 
\v{E) u*{-E)) 


(D9) 









9 


that solves the homogeneous equation 

(^-tla,-^+E + mWW'^^X{E)=Q. (DIO) 

The functions u and v are determined by 

7ou" + (e^ + — *)u = 0, (v = ieu — ^ou', (Dll) 


(D12) 

(D13) 


(DM) 


e = {E + iTT^)/jo, C = ( 27 r/ 7 o)(a|/ 3 )*, 

= ICI = i(’’'/7o)(7i - 72 ). 

The solution is 
u{E) = 

Cv{E) = [/(i _ 1 i. 

Finally, the scattering matrix has the dyadic form 

Sn^{E,E')\E>E' =-^n(E)^l,{-E'), (D15) 

^{E) = (27r/7o)e-^'/8i^t ME)\ 


independent of the rapidity 70 of the fermion-parity 
switch. 

As a check, we can directly compute the transferred 
charge in the adiabatic limit from Eq. (301. Substitution 


of the frozen scattering matrix at the Fermi level, 

So = I + 2TriW\Eoa^ - inWW^y^W, (D19) 
gives the charge 


Qadiabatic — . [ dEQT^lSQCTz 

47r ^ dEo 


- / dEo Tr {Eoa^ + iTvWW^^y^Wa^W^^ 
' J — 00 

• [Eoa, - iTvWW^y^a, 

/ OO 

dEoiEl+^y^)-^ 

-00 


= eTZ. 


(D20) 


3. Transferred charge 


4. Relation of the reduction factor to the Andreev 
conductance 


Because the scattering matrix is still of rank-one, a sin¬ 
gle quasiparticle is transferred as a result of the fermion- 
parity switch, irrespective of the number of channels N 
in the metal probe. The charge expectation value of this 
quasiparticle is given by 

Q = e/ —r{E)ay{E) 


/o 


2 !^ E 

7o 


/o 

^e-^^/4((a|a)-(/3|/3)) 
7o 




\viE) 


pOO 

X dE{\u{Ey-\v{Ey). (D17) 
Jo 


Comparison with Eq. (28) shows that the transferred 


charge for a multi-channel contact differs from that in 
the single-channel case by a reduction factor 


7^ = 


(ala) - (/3|/3) 


(a|a) - (/3|/3) 


C((«|a) + (/3|/3))2-4|(a|/3)P 


e[o,i], (DI 8 ) 


The charge reduction factor TZ from Eq. (D18) is a 


property of the coupling matrix of the normal-metal 
probe to the Josephson junction. It can be expressed 
in terms of an independently measurable quantity, the 
Andreev conductance. 

When the normal-metal probe is biased at a voltage V, 
a current I is driven into the grounded superconductor by 
the process of Andreev reflection. The Andreev conduc¬ 
tance Ga = limv_>o dljdV is related to the scattering 
matrix So at the Fermi level by 


Ga = ^ Tr (1 - Soa^S^a^. (D21) 


Near the level crossing a resonant peak appears in Ga 
as a function of Eq, with the Lorentzian line shape 


4e^ 

Ga = -^ „2 ^ 2 - 2 ( 1 -^ )■ 

h Eq + ^ ' 


(D22) 


The resonant peak height of (4e^//i)(l — TZ^) directly de¬ 
termines the charge reduction factor. 
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